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We determine the second order endpoint of the line of first order phase transitions, which occur in the light 
quark mass regime of 3-flavour QCD at finite temperature, and analyze universal properties of this chiral critical 
point. A detailed analysis of Binder cumulants and the joint probability distributions of energy like and ordering- 
field like observables confirms that the chiral critical point belongs to the universality class of the 3d Ising model. 
From a calculation with improved gauge and staggered fermion actions we estimate that the transition is first 
order for pseudo-scalar meson masses less than about 200 MeV. 



1. INTRODUCTION 

The phase transition in 3-flavour QCD is first 
order in the chiral limit (m = 0). This first order 
transition will continue to persist for small but 
non-zero values of to up to a critical value, m, 
of the quark mass. It is expected that QCD at 
this chiral critical point belongs to the universal- 
ity class of the 3d Ising model We will present 
here evidence for this scenario 

The lattice formulation of QCD with three 
degenerate quark masses depends on the bare 
quark mass m and the gauge coupling (3 = 
The thermodynamics is described in terms of the 
partition function Z{f3,m) = / 2?[/e-S({;7},/3,m)_ 
The Euclidean action, 5, is given in terms of the 
pure gauge (Sg) and the fermion action (Sp) 

S{{U}, f3, to) = PSg{{U}) - Sf{{U}, to) , (1) 

where Sp = 0.75TrlnM is expressed in terms 
of the fermion matrix M{{U}). We will mainly 
use standard staggered fermions and the standard 
Wilson gauge actions. 

The universal properties at the chiral critical 
point are controlled by an effective Hamiltonian, 



(2) 



with a global Z{2) symmetry. It controls the criti- 
cal behaviour at {(Scifh), rfi) and can be expressed 
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in terms of two operators f , A^, i.e. the energy- 
like and ordering-field like operators that couple 
to two relevant scaling fields r and ^. As the sym- 
metry of Tief / is not shared in any obvious way by 
the QCD Lagrangian one may expect that in the 
vicinity of the chiral critical point the operators 
and couplings appearing in the QCD Lagrangian 
are linear combinations of £, Ai and r, ^ respec- 
tively. One may use a linear ansatz 

r = /3 — /3c + ^ (to — m) , 

^ = TO-TO + B(/3-/3e), (3) 

as well as 

£ ^ Sg + r ^pi; , M = ^ptp + s Sg ■ (4) 

Here = 0.75 TrAf^^ denotes the chiral con- 
densate expressed in terms of bosonic observ- 
ables. As Sg, V'V' or related observables like the 
Polyakov loop L are mixtures of £ and A4, the 
corresponding susceptibilities will all receive con- 
tributions from fiuctuations of £ as well as A4. 
Asymptotically therefore all of them will show 
identical finite size scaling behaviour which will 
be dominated by the largest universal exponent 
7/1/. Unfortunately this is quite similar for the 
3d-Z{2) and 0{N) spin models, so that a finite 
size scaling analysis of susceptibilities is not a 
good indicator for the universality class. The sit- 
uation is different for cumulants constructed from 
linear combinations of tpip and Sg, 



Bi{x) 



{{SM{x)y 



M{x) =tp^lj + xSG. (5) 
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Here SX denotes X — {X). For all values of x 
different from 1/r the cumulants behave asymp- 
totically like the Binder cumulant for the order 
parameter; cumulants calculated on different size 
lattices for different quark masses will intersect at 
some value of the quark mass. In the infinite vol- 
ume limit these intersection points will converge 
to a universal value which is characteristic for the 
universality class, in fact, it is quite different for 
the classes of 3d-Z{2) and 0{N) symmetric spin 
models; e.g. = 1.604 for Z(2) [|, 1.242(2) for 
0(2) g and 1.092(3) for 0(4) @. The cumulants 
B^^x) thus seem to be appropriate observables to 
locate the chiral critical point as well as to de- 
termine its universality class without knowing in 
detail the correct scaling fields. 

2. LOCATING THE CRITICAL POINT 

We make use of the finite size scaling prop- 
erties of Binder cumulants Bi{x) evaluated at 
f3pc{m). From previous studies one knows that 
the endpoint in 3-flavour QCD is located close to 
TO = 0.035 We have performed calculations 
on lattices of size x 4, with = 8, 12 and 
16. We have used four values of the quark mass 
in the interval to G [0.03, 0.04] and for each of 
these masses we calculated thermodynamic ob- 
servables for 3 to 4 different values of the gauge 
coupling (3. In general we collected for each pair 
of couplings (1 — 3) • 10^ configurations generated 
with the hybrid-R algorithm^. Interpolations be- 
tween results from different /^-values have been 
performed using the Ferrenberg-Swendsen multi- 
histogram technique. The cumulant of the chi- 
ral condensate, -84(0), is shown in Fig. |l]. We 
note that the cumulants calculated on different 
size lattices intersect at a quark mass close to 
TO — 0.035. The value of ^4(0) at the intersection 
point is compatible with the universal value of the 
Binder cumulant for the 3d Ising model. We have 
fitted the cumulants on a given lattice size using 
a linear ansatz in the quark mass. From this we 
find for the critical mass and the cumulant, 



TO = 0.033(1), ^4(0) = 1.639(24) , 



(6) 
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Figure 1. The cumulant -84(0) defined in Eq. 5. 

This result strongly suggests that the chiral crit- 
ical point indeed belongs to the universality class 
of the 3d Ising model. As discussed above the in- 
tersection point of Binder cumulants constructed 
from M{x) according to Eq. ^ will be independent 
of the choice of x in the infinite volume limit. The 
optimal choice a; = s, will minimize the finite vol- 
ume effects, which are quadratic in A = a; — s. 
The determination of the a;-value for which the 
intersection point of the Binder cumulant is clos- 
est to the Z{2) value provides a finite volume es- 
timate for the mixing parameter s. One finds 
Smin = 0.430(23). 

3. THE MIXING PARAMETERS 

The parameters r and s are fixed by demanding 
that A4 should obey basic properties of the order 
parameter for spontaneous symmetry breaking at 
the critical point. We obtain two conditions for 
the mixing parameters r, s in terms of the pa- 
rameter B and expectation values of Sq, 4''4'-i 



= -B. 



{{5SgY) - B{SHj SSg) 



(7) 



^We used trajectories of length r = 0.675 generated with 
a discrete step size St = 0.015. 



The parameter B is needed to define lines of con- 
stant ^. The line of first order phase transitions 
defines the zero external field line (^ = 0) of the 
effective Hamiltonian. We thus can extract B 
from the quark mass dependence of the pseudo- 
critical couplings. Using Eq. |^ one obtains 

/ d/3c(TO 



B 



\ dm 



(8) 
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Figure 2. Joint probability distribution of the 
fluctuations in the ordering-field hke and energy 
Uke operators, which have been constructed ac- 
cording to Eq. 1^ with r = 0.55 and s = 0.43. 

Knowing B we also know r and can construct s 
using Eq. ^. A first estimate may be given using 
our data on the largest lattice (f6^x4). The slope 
of Ppdm) defines r~^. We estimate r = 0.5f (2) 
from a straight line fit. The mixing parameter 
r is large and definitely non-zero. The impor- 
tance of choosing the correct mixing parameter r 
becomes apparent from an analysis of joint prob- 
ability distributions for S£ and SA4. These are 
little affected by changes of s but they strongly 
depend on r. In order to improve our determi- 
nation of r, we demand that fiuctuations in Ai 
vanish for any fixed value of S£. This maximizes 
the Z{2) symmetry of the contour plot shown in 
Fig.||. One obtains in this way r = 0.550(7). A 
determination of the parameter s from Eq. |^ then 
yields s — 0.41(51). Although this value has large 
errors it is consistent with the result found from 
the intersection point of the Binder cumulant. As 
a best estimate we therefore obtain 



r = 0.550 ±0.007, s = 0.430 ± 0.023 . 



4. THE PHYSICAL SCALE 



(9) 



In order to determine a physical scale for the 
chiral critical point we have calculated hadron 
masses on a IG'^ lattice at (/3c, m). For the pseudo- 
scalar meson mass we find nips/Tc = 1.853(1). 
Using estimates for the critical temperature we 
thus estimate for the pseudo-scalar meson mass 
at the chiral critical point nips — 290 MeV. The 



entire analysis of the critical point discussed so 
far has been performed with unimproved gauge 
and fermion actions on rather coarse lattices. Im- 
proved actions are not expected to modify the 
results on the universal properties of the criti- 
cal point. They may, however, well influence the 
quantitative determination of the critical point. 
We therefore have investigated the critical point 
also in calculations with improved gauge and 
staggered fermion actions, the p4-action 
We have performed calculations on a lattice of 
size 12^ X 4 with a bare quark mass m = 0.005 
and on a lattice of size 16'^ x 4 with m = 0.01. 
Again a Ferrenberg-Swendsen reweighting is used 
to determine the pseudo-critical couplings and 
Binder cumulants at these couplings. The smaller 
quark mass leads to a Binder cumulant below 
the Z{2) value (^4(0) = 1.31(12)) and for the 
larger quark mass, the Binder cumulant is signif- 
icantly above the Z{2) value (^4(0) = 2.14(10)). 
We thus have obtained an upper and lower limit 
for the critical point, suggesting a critical bare 
quark mass of ?fi = 0.0075(25). Extrapolat- 
ing the meson masses calculated in to this 
value of the bare quark mass, we estimate for 
the pseudo-scalar meson mass at the critical point 
TTips ~ 192(25) MeV. The critical mass thus is 
considerably smaller than estimated from calcu- 
lations with unimproved gauge and fermion ac- 
tions. 
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